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1. Introduction 

Photons emitted by the quark-gluon plasma (QGP) created at the Relativistic Heavy Ion 
ColUder (RHIC), as well as in future heavy ion experiments at the Large Hadron Collider 
(LHC), are expected to be an important probe of the plasma physics. Because of the finite 
extent of the plasma ball and of the weakness of the electromagnetic coupling constant, once 
a photon is produced it escapes essentially undisturbed, thus carrying valuable information 
about local properties of the plasma at the point of emission [1]. 

At weak coupling, photon production in Quantum Chromodynamics (QCD) can be stud- 
ied by means of perturbation theory [2]. However, results from RHIC seemingly indicate that 
the plasma created there does not behave as an almost free gas of quarks and gluons, but 
rather as a strongly coupled liquid [4] . Consequently, understanding the spectrum and rate of 
photon emission requires a non-perturbative calculation. This is currently beyond the scope 
of analytical methods, and performing this calculation on the lattice is problematic due to 
the real-time, Lorentzian-signature nature of the physics involved. 

In view of the above, it is interesting to investigate photon production in the context of the 
gauge/gravity correspondence [5]. Although the gravity dual of QCD itself is not known, the 
duals of a large class of finitc-tcmperature, strongly coupled, large- A'c SU{Nc) gauge theories 
are known. One may hope that generic, universal features exhibited by such theories may be 
useful indications as to the behaviour of QCD, at least in the regime of strong coupling. With 
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this in mind, ref. [6] studied photon production at strong couphng in = 4 super Yang-Mills 
(SYM), a four-dimensional theory with massless matter in the adjoint representation, and 
ref. [7] calculated the rate of photon emission in a four-dimensional theory with massless 
quarks. One purpose of this paper is to study the effect of a non-zero quark mass on photon 
production. 

In perturbative, weak-coupling calculations it is usually assumed that the temperature T 
is high enough so that the thermal quark mass, of order \/AT (with A = the 't Hooft 

coupling) is much larger than the bare, zero-temperature quark mass Mq, and so the latter is 
effectively set to zero. Regardless of the precise temperature range in which this is applicable, 
the very nature of this approximation assumes that quarks are well defined quasi-particles to 
which physical properties such as a thermal mass can be meaningfully assigned. This is true 
at sufficiently weak values of the coupling, but it is not necessarily true at strong coupling. 
Well defined quasi-particle excitations certainly do not seem to exist in the phases of interest 
of the holographic theories that we will consider [8, 9], and they might not exist either in a 
QCD plasma at temperatures just above deconfinement. Under these circumstances, the bare 
quark mass is better thought of as a microscopic parameter on which the physics depends, 
but without a direct interpretation as the mass of a physical quasi-particle. Examples of 
physical quantities that depend on the quark mass through the dimensionless combination 
MJT include the entropy density [8, 9], the shear viscosity [10], etc. We will see that the 
same is true for observables related to photon production. 

We will consider finitc-tcmperature, SU(N^) super Yang-Mills (SYM) theories coupled 
to Nf flavours of fundamental matter whose bare mass is an adjustable parameter. The 
fundamental matter includes both fermions and scalars, which we will refer to collectively 
as 'quarks'. Our goal is to study photon production as a function of the quark mass and 
the photon frequency, to leading order in the electromagnetic coupling constant but non- 
perturbatively in the SU{N^) interactions. Note that, in addition to the quarks, the matter 
content in these theories also includes adjoint fields. Since we would like to model as closely as 
possible QCD, which has no adjoint matter, we will assume that the photon couples directly 
only to the quarks, i.e., wc will assign vanishing electric charge to the adjoint matter. 

In the limit of a small number of flavours, A^f <C A';., the theories above possess a simple 
dual description in terms of Nf Dq-brane probes in the gravitational background of Dp- 
branes [11]. At temperatures high enough so that the gluons (and the adjoint matter) are 
deconfined, the gravitational background contains a black hole [12], represented by the shaded 
blob in fig. 1. In this deconfined phase, the fundamental matter may still be in two different 
phases separated by a first-order phase transition [9].^ 

Prom the viewpoint of the holographic description, the basic physics behind this transition 
is easily understood. The asymptotic distance between the Dq-branes and the black hole is 
proportional to the quark mass, whereas the size of the black hole horizon is proportional 
to the temperature. Thus, for sufficiently small T/M^, the Dq-branes are deformed by the 

'^Specific examples of this transition were originally seen in [13, 14]. Aspects of these transitions were 
independently studied in the D3/D7 system in [15, 16], and in a slightly different framework in [17]. 
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Figure 1: Possible embeddings for Dq-brane probes in the background of black Dp-branes. The 
critical embedding is skipped over by the phase transition. 

gravitational attraction of the black hole, but remain entirely outside the horizon in what 
we call a 'Minkowski embedding' (see fig. 1). However, above a critical temperature Tf^^, 
the gravitational force overcomes the tension of the branes and these are pulled into the 
horzion. We refer to such configurations as 'black hole embeddings'. In between these two 
types of embeddings there exists a 'critical embedding' in which the branes just 'touch the 
horizon at a point'. However, thermodynamic considerations reveal that a first-order phase 
transition occurs between a Minkowski and a black hole embedding. In other words, the 
critical embedding is skipped over by the phase transition, and near-critical embeddings turn 
out to be metastable or unstable. 

In the dual field theory, this phase transition is exemplified by discontinuities in physical 
quantities such as, for example, the quark condensate or the contribution of the fundamental 
matter to the entropy density. However, the most striking feature of this phase transition 
is found in the spectrum of physical excitations of the fundamental matter. In the low- 
temperature, Minkowski phase the spectrum is gapped and contains a discrete set of deeply 
bound mesons {i.e., quark-antiquark bound states) with masses of order M^^^ ~ M^/y/X. 
These mesons are dual to excitations supported on the probe branes - see, e.g., [18, 19, 20]. 
In addition, the Minkowski-phase spectrum also contains well defined, quark-like excitations 
described by strings stretching between the tip of the branes and the horizon. These have 
masses of order and are therefore parametrically heavier than the mesons. Both sets of 
excitations are absolutely stable in the large- A*";., strong coupling limit under consideration. 

In the high-temperature, black hole phase stable mesons cease to exist. Rather one finds 
a continuous and gapless spectrum of excitations [8, 9, 21, 22]. Hence the mesons dissociate or 
'melt' at the first order phase transition at Tf^^, and beyond it no well defined, quasi-particle 
notion of an individual quark exists, since a string stretching between any point on the branes 
and the horizon will quickly fall through the horizon. In the gauge theory this corresponds 
to the fact that any localised quark charge will quickly spread across the entire plasma, thus 
loosing its identity. 

In this paper we will study photon production in the black hole phase. We will see that 



- 3 - 



this is proportional to NfN^, i.e., to the number of electricahy charged degrees of freedom, as 
one would expect. For = our result coincides with that of [6], except for the fact that in 
that reference the overall normalisation is set by N^, reflecting the fact that electric charge is 
carried by adjoint degrees of freedom. We stress that the photon production we will calculate 
is not a correction to an 0(A^_f )-rcsult but the leading result in the large- A^^ expansion, since 
in our case the adjoint matter is electrically neutral. We will see that the gravity picture 
provides a simple, geometric explanation for the agreement (up to normalisation) between 
our Mq = result and that of [6]. 

We will not explicitly study photon production in the Minkowski phase, since we have a 
clear expectation of what the result would be. Indeed, as explained above, in this phase the 
spectrum consists of a discrete set of stable modes, so we expect the spectral function (2.4) 
to be given by a sum of delta functions supported at appropriate values of the momentum. 
We will confirm this expectation indirectly by calculating the spectral function for black hole 
embeddings near the critical embedding (see fig. 1). 

In order to study photon production, one must in principle introduce a dynamical photon 
that couples appropriately to the SU{N^) fundamental matter and construct the dual of the 
SU{N,) X U{1) EM gauge theory thus obtained. However, we will see that, to leading order 
in the electromagnetic coupling constant, photon production in the SU{N^) x U{1)e,m theory 
is completely determined by correlators in the original SU (iV^) theory with no dynamical 
photon. This observation allows the holographic calculation of the desired correlators. 

This paper is organised as follows. In sec. 2 we discuss how to calculate the rate of 
photon production in a gauge theory at finite temperature. In sec. 3 we explain how to use 
the holographic description in terms of a Dp/Dq system to study photon production in the 
dual gauge theory. In sees. 4 and 5 we present our results for the D3/D7 and the D4/D6 
systems, respectively. We conclude in sec. 6 with a discussion of our results and their possible 
implications for heavy ion collision experiments. 

2. Photon emission in a hot gauge theory 

Consider a (d+l)-dimensional SU{N^) gauge theory with Nf flavours of fundamental fermions 

and scalars a = 1,...,N{, to which we will collectively refer as 'quarks'. Imagine 
coupling this theory to electromagnetism by introducing and additional, dynamical, Abelian 
gauge field, the photon, that couples to the fundamental fields with strength (electric charge) 
e.^ This new theory may be constructed by adding a kinetic term for the photon and replacing 
the S'C/(Ai'(.)-covariant derivative by = — ieA^ when acting on the fundamental 
fields, with n = 0, . . . ,d. In this way we obtain an SU{N^) x U{1)-em gauge theory with 

^For simplicity we assume that all the fundamental fields have equal mass and carry equal electric charge, 
but our results can be generalised straightforwardly. Later we will consider SU{Nc) theories that also contain 
adjoint matter. Since our goal is to model QCD, which has no such matter, we will assign vanishing electric 
charge to the adjoint fields. 
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(a) (b) 

Figure 2: Diagrams contributing to the two-point function of electromagnetic currents (2.5). The 
external line corresponds to a photon of momentum k. As explained in the text, to leading order in 
the electromagnetic coupling constant only SU{Nc) fields 'run' in the loops represented by the shaded 
blobs. 



Lagrangean 

C = CsuiN.)-l^'. + eA^j;'', (2.1) 
where J^^y = d/^Av — di,Aix and the electromagnetic current is given by 

J^- = + (P^$) - iV^^r $ . (2.2) 

A sum over flavour and colour indices is implicit in this formula. Note that the photon field 
enters this current through the covariant derivative acting on the scalars. 

In thermal equilibrium, the differential photon emission rate per unit time and volume, 
at leading order in the electromagnetic coupling constant e, is then given by [23] 

^ = J2^\''^^^"^ i:i)(k)4.)(k) X..(/c)|,o=, , (2.3) 

V / II g—i 

where k = (fc^jk) is the photon null momentum, 

X^.(fc) = -2ImG^,(fc) (2.4) 

is the spectral density, and 

GlAk) = -iJ d'+'xe-"^-^e{x'){[j;^ix),Jrm) (2-5) 

is the retarded correlator of two electromagnetic currents, whose diagramatic representation 
(including photon fields as external legs) is given in fig. 2. Finally, 

n.{k^) = (2.6) 
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is the standard Bose-Einstein distribution function. Without loss of generaUty we will assume 
that k points in the x^-direction, and we will denote by x*, z = 2, . . . , d the remaining spatial 
directions. The polarisation vectors e^^^ may be chosen to be unit spatial vectors orthogonal 
to k. Transversality of the correlator implies that the sum over polarisation vectors in (2.3) 
may be replaced by the trace of the spectral function, i.e., 

d-l 

X^,{k) = r^^^xAk) = Y.is)i^yis){^)xA^) , (2-7) 

s=l 

and rotational invariance ensures that this is given by 

x';(A;) = -2(cZ-l)ImG^(A;), (2.8) 

where 

G''{k)^j^5^^Gfj{k). (2.9) 
The trace of the spectral function also determines the electric conductivity as 

„2 1 



a 



— lim x^aik)\^o u ■ (2-10) 

Thus in order to study photon production we must in principle calculate the two-point 
function (2.5) in the SU{N^) x C/(1)bm theory. However, as noted in [6], to leading order in 
the electromagnetic coupling constant this reduces to a calculation purely within the original 
SU (A^c) theory; this is the key observation that will allow us to perform this calculation using 
the dual gravitational description, since the dual of the SU^N^) x U (1)em theory is unknown. 
To see this, note first that the terms in the electromagnetic current (2.2) proportional to the 
photon field (implicit in the covariant derivative) lead to higher-order contributions in e to 
the correlator (2.5), and can thus be ignored to leading order in e. Second, observe that the 
two-point function of the remaining terms in the current can be calculated in the SU{N^) 
theory, since again the effects of including the dynamical photon are of higher order in e. 
Diagramatically, this means that no photon fields are present in the shaded blobs in fig. 2. 
An additional observation is the fact that the tadpole diagram in fig. 2 has no imaginary part 
and hence does not contribute to the desired spectral function. 

We therefore conclude that, to leading order in the electromagnetic coupling constant, 
photon production in an SU{N^) x U{1)em theory is completely determined by the two-point 
function of the electromagnetic current in the SU (N^) theory. In the rest of the paper we will 
calculate this correlator in SU{N^) SYM theories coupled to fundamental matter. 

3. Holographic description 

In the decoupling limit, the black Dp-brane solution in the string frame takes the form 



e* = i?^, Coi...p = H-', (3.1) 



3-p 
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where x = (x^, . . . 



fir) = 1 




(3.2) 



and L is a length scale (the AdS radius in the case p = 3). The horizon lies at r = tq. As 
usual, regularity of the Euclidean section, obtained through — > ix^, requires that x% be 
identified with period 



In this paper it will be useful to work with a different radial coordinate u related to r through 



in terms of which / = 1 — 4tt^ and the horizon lies at n = 1/2. 

According to the gauge/gravity correspondence, string theory on the background above is 
dual to {p+ l)-dimensional, maximally supersymmetric Yang- Mills theory at temperature T. 
In some cases one periodically identifies some of the 'Poincare' directions x in order to render 
the theory effectively lower-dimensional at low energies; a prototypical example is that of a D4- 
brane with one compact space direction. Under these circumstances a different background 
with no black hole may describe the low-temperature physics, and a phase transition may 
occur as T increases [12]. In the gauge theory this is typically a confinement/deconfinement 
phase transition for the gluonic (and adjoint) degrees of freedom. Throughout this paper we 
assume that T is high enough, in which case the gauge theory is in the deconfined phase and 
the appropriate gravitational background is (3.1). 

Consider now Nf <^ coincident Dq-brane probes that share d spacelike Poincare di- 
rections with the background Dp-branes and wrap an iS" inside the S^~^. We will assume 
that the Dq-branes also extend along the radial direction, so that q = d + n + 1. We will 
denote by x^, with n = 0, . . . ,d, the directions common to both branes. In the gauge theory 
the Dq-branes correspond to introducing Nf flavours of fundamental matter of equal mass 
that propagate along a {d + l)-dimensional defect. To ensure stability, we will assume that 
the Dp/Dq intersection under consideration is supersymmetric at zero temperature. Under 
these conditions the Ramond-Ramond field sourced by the Dp-branes does not couple to the 
Dq-branes. Two cases of special interest here are the D3/D7 (n = 3) [13] and the D4/D6 
{n = 2) [14] systems. If one of the D4 directions is compact, then both cases can effectively be 
thought of as describing the dynamics of a four-dimensional gauge theory with fundamental 
matter. 

The U{Nf) gauge symmetry on the Dq-branes is a global, flavour symmetry of the gauge 
theory. If this theory is coupled to electromagnetism as explained in the previous section 
(in particular, if all flavours are assigned equal electric charge e), then the electromagnetic 
current J^"^' is the conserved current associated to the overall U{1) C U{N(). In order to 




(3.3) 




1 



7-p 



(3.4) 
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study photon emission to leading order in e it suffices to evaluate the correlator (2.5) in the 
SU{N^) gauge theory with no dynamical photon. 

At strong 't Hooft coupling and large N^, this correlator is easily calculated holographi- 
cally. Indeed, global symmetries of the gauge theory are in one-to-one correspondence with 
gauge symmetries on the gravity side, and each conserved current of the gauge theory is dual 
to a gauge field on the gravity side [24]. Let A^, m = 0, . . . ,q, be the guage field associated 
to the overall C/(l) gauge symmetry on the Dq-branes. Upon dimensional reduction on the 
n-sphere wrapped by the Dq-branes, gives rise to a massless gauge field {A^, A^}, n mass- 
less scalars, and a tower of massive Kaluza-Klein (KK) modes. All these fields propagate on 
the d + 2 non-compact dimensions of the Dq-branes. We will work in the gauge A^ = 0, and 
we will consistently set to zero the scalars and the higher KK modes, since they are not of 
interest here. The gauge field Ai^ is the desired dual to the conserved electromagnetic current 
J™ of the gauge theory; just like the quarks axe localised on a defect inside the gauge theory, 
the dual gauge field is localised on the Dq-branes inside the bulk spacetime. 

According to the gauge/gravity correspondence, correlation functions of J™ can be calcu- 
lated by varying the string partition function with respect to the value of A^ at the boundary 
of the spacetime (3.1). Under the present circumstances, the string partition function reduces 
to e*'^, where S is the sum of the supergravity action and the effective action for the Dq-branes. 
Since does not enter the supergravity action, the form of this action will not be needed 
here. Moreover, the Wess-Zumino part of the Dq-branes action does not contribute for the 
brane orientations and the gauge field polarisations considered in this paper. Therefore, it 
suffices to consider the Dirac-Born-Infeld part of the Dq-branes action: 

S = -NfTj,^[ (P^^x e-^^- det(5 + ^-kIIF) , (3.5) 

JDq 

where g is the induced metric on the Dq-branes, F = dA is the overall U{1) field strength, 
Toq = 1/ {2TTisYgs^s is the Dq-brane tension, and ig and gg are the string length and coupling 
constant, respectively. Since we are only interested in the two-point function (2.5), we only 
need to keep terms up to quadratic order in the gauge field. We will therefore work with the 
simpler action 

S = -N,T^, d'l+'x e-'^y/^d^ (^1 + ^^F^^ , (3.6) 
where = FmnF"'"'. 

The Dq-branes wrap an n-sphere in the directions transverse to the Dp-branes, so it is 
convenient to write the metric on S^~p in adapted coordinates as 

dnl_p = d9^ + sin^ 9 d^l + cos^ 9 dQ\_p_^ . (3.7) 

Setting il) = cos 9, the Dq-branes embedding may then be specified as ■0 = i^iu). Asymptoti- 
cally, ie., as r — oo or u — 0, this behaves as [8, 9] 

2 2n 

'4>{u) ~ mu''-p + cu''-p + • • • , (3-8) 
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where m and c are proportional to the quark mass and condensate, respectively (see below). 

Since the gauge field enters the action (3.6) quadratically, turning on A induces a cor- 
rection Sip of order on the branes embedding. When the action is evaluated on-shell, this 
generates a boundary-term contribution of the form 



(3.9) 

boundary 



where £. is the Dq-branes Lagrangean density and ijj' = dip/du. Since Sip = 0{A^), this has in 
principle the right form to contribute to the two-point function (2.5). However, as explained 
in [9], this contribution is proportional to 6m, and hence vanishes because the quark mass 
must be kept fixed in computing the correlator. 

We thus conclude that, in order to calculate the photon emission rate, we may consistently 
proceed by first determining the Dq-branes embedding in the absence of the gauge field, and 
then solving for the gauge field on that embedding. As explained above, we will set to zero 
the components of the gauge field on the n-sphere wrapped by the Dq-branes. Moreover, 
following [25], we will choose the gauge = 0, Fourier-decompose the d + 1 remaining 
components of the gauge field as 

A,{x', ^,u)=J e-^'=°-°+^>^- A,{k', k, u) , (3.10) 

and choose k to point in the a;^-direction. Under these circumstances the equations of motion 
for the Aj-components of the gauge field (i = 2, . . . , d) that follow from the action (3.6) 
decouple from each other and from those for Aq,Ai. Setting Ai = A(k,u), the equation of 
motion takes the form 



5, 



Qiu)duA +k1 l-f{u) P{u)A = 0, (3.11) 



where we have made use of the fact that |k| = The functions P{u),Q{u) depend on the 
embedding ip{u) and will be specified below. Eq. (3.11) may be consistently derived from the 
'reduced' action 

S = -^o^J dx^d'^^du [-P{u){dQAf + fP{u){diAf + Q{u){duAf] , (3.12) 

obtained from (3.6) by integrating over the n-sphere and setting diA = 0. Wc will sec 
below that the normalisation constant in front of this action scales as AAdq ^ NfN^T'^~^ and 
determines the overall magnitude of photon production. 

According to the prescription in [26, 27], the retarded correlator (2.9) is then given by 

G-{u) = -AT,, hm ^Qi^)A*{u. u)d,^A{u,u) 

^ ' A*{uj,d)A{io,Q) ' ^ ' 

where u = /2TiT and A{uj,u) is a solution of the equation of motion (3.11) obeying the 
incoming- wave boundary condition at the horizon. The imaginary part of G^(u;), which is all 
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we need, is tt-independent, and it is convenient to calculate it at the horizon instead of at the 
boundary: 

ImG^(u;) = "^"^ Im lim 2Q{u)A*{oj,u)duA{u},u) . (3.14) 

|A(a;,0)| u^^ 

Note that, although not explicitly indicated, several quantities above depend not only on the 
photon frequency u but also on the quark mass, since the equation of motion (3.11) that 
determines A depends on the branes embedding through the functions P and Q. 



4. The D3/D7 system 

Here we will specialise the above discussion to the D3/D7 system. This intersection is sum- 
marised by the array 

0123456789 
D3: X X X X (4.1) 
D7: xxxxxxxx 

Of course, this is an interesting system because d = 3, i.e., both the gluons and the funda- 
mental fields in the gauge theory propagate in 3 + 1 dimensions. 
The metric (3.1) for black D3-branes takes the form 

ds^ = ^ {-fdxl + dx2) + + L'dnl , (4.2) 



where 



4 



L"" = iirgsNJt , / = 1 - , ro = ttTL' . (4.3) 



r 



In this case the dimensionless coordinate (3.4) is^ u = rQ/2r^, in terms of which the metric 
becomes 

ds' = {-fdxl + + + ^"dnl . (4.4) 

Specifying the D7-branes embedding through = ip{u), as explained above, the induced 
metric takes the form 

n (7rTL)2 „ L2 (1 - + 4n2/V''^) 9 9. 9. 9 . . 

dslr = {-fdxl + d^') + 4^2^(1 _ ^2) ' du' + ^'(1 " " (4-5) 

In this case the asymptotic behaviour of -0 is^ 

i}{u)=m v}!"^ + c u^l'^ + • • • , (4.6) 



'This is related to the coordinate used in [25] as 2Mhcre = M[25] • 

*In this hmit, the coordinate u used here and the coordinate p used in [8, 9] are related through = 1/t 
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where the dimensionless constants m and c are related to the quark mass and condensate, 
respectively, through [8, 9]: 

(O) = -2^/^7r^£lNtTj,r4c= -lVXN(N,T^c, (4.8) 

8 

with A = 5ym-^c = "^T^gs^c the 't Hooft coupling. The operator O is a supersymmetric version 
of the quark bilinear, and it takes the schematic form 

O = + + Mq^'f^ , (4.9) 

where X is one of the adjoint scalars. We will loosely refer to its expectation value as the 
'quark condensate'. A detailed discussion of this operator, including a precise definition, can 
be found in the appendix of ref. [28]. 

Eq. (4.7) implies the relation m = M/T between the dimensionless quantity m, the 
temperature T and the mass scale 

7rL2 27r ' 

where M^es is the mass of the lightest meson (i.e., the mass gap) in the discrete meson 
spectrum at zero temperature [18, 19, 20, 14]. M determines the scale of the temperature of 
the phase transition for the fundamental degrees of freedom, Tf^n ~ M, since the latter takes 
place at m ~ 1. 

The functions P{u), Q{u) occurring in the reduced action (3.12) are given in this case by 



P(u) = (^- ^^^)V^-^'' + ^>i f^ Qiu) = , - (4 11) 

^^"^ u/(27rr)2 ' 2Vl-V'2+4uW ^^■''> 



and the normalisation constant is 



where 

AT 32 



AT,, = !^^J^J^ = IxN.N^T^ (4.13) 



is the normalisation constant introduced in [8, 9] and il^ = 27r^ is the volume of a unit 
three-sphere. 

Analysis of eq. (3.11) near u = 1/2 reveals that the solution obeying the incoming-wave 
boundary condition behaves as A ~ (1 — 2n)~"^/^ near the horizon. We will thus seek a 
solution of the form 

A{u,u) = {1 - 2u)-''^/^{l + 2u)-^/^ F{u,u) , (4.14) 
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where F is a regular function of u and the second factor has been exphcitly extracted for 
convenience. Since the differential equation solved by A is linear, its overall normalisation 
is arbitrary, and it cancels out in eqs. (3.13), (3.14). We choose to fix it by setting F equal 
to unity at the horizon, i.e., F{uj, ^) = 1. Using the form of the solution (4.14) and this 
normalisation choice for F, the limit on the right-hand side of (3.14) is readily calculated 
with the result 

lim 2Q{u)A*{u,u)duA{u,u) = ^ , (4.15) 

where ipo = V'(l/2) is the value of at the horizon. Thus the trace of the spectral density, 
eq. (2.8), takes the form 

X';(c.) = -4ImG-(^) = 8^-0, (1 - i^lf J' . (4.16) 
4.1 Massless quarks 

Massless quarks correspond to the equatorial embedding ^ = 0. In this case 

and the equation of motion (3.11) becomes 

(1 + 2uf{l - 2uf A" - 8n(l + 2u){l - 2u) A' + %iJu A = Q, (4.18) 

where uj = /2'kT and a prime denotes differentiation with respect to u. Upon changing 
variables as described in footnote 3, the equation of motion above agrees precisely with 
eq. (5.5d) in [25]. The reason for this is that, for massless quarks, the induced metric (4.5) 
on the D7-branes is exactly AdS^ x S^' . After reduction on the , the quadratic action for 
the gauge field on the D7-brancs is then, up to an overall normalisation, exactly the same as 
that considered in [25] for a gauge field in AdS^. The normalisation constant we obtain scales 
as TVAd? ~ -^f-^c) whereas that of [25] scales as A^^^. This merely reflects the difference in the 
number of electrically charged degrees of freedom in the corresponding dual gauge theories. 

The solution of (4.18) satisfying the incoming- wave boundary condition was found in 
ref. [6]: 

A{u,u) = (l-2«)-^-/2(l + 2n)-/SFi (^a,6,c;i^^ , (4.19) 

where 

a = l-{l + i)'^, b=-{l + i)'^, c=l-iu}, (4.20) 
and 2-Pi is a hyper geometric function. Substituting into eq. (4.16) we obtain the result 

= , Jf^-;" (4.21) 
2-|2Fi(a,6,c;i)| 

for the trace of the spectral density. A plot of the dimensionless combination x^ti{'^) I ^^tit ^ 
is given by the top, solid, red curve in fig. 3. 
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Figure 3: D3/D7 system: Trace of the spectral function as a function of lo for (from top to 
bottom on the left-hand side) m = {0,0.6,0.85,0.93,1.15,1.25,1.306}, or equivalently for t/jq = 
{0,0.37,0.53,0.58,0.75,0.85,0.941}. The last value corresponds to that at which the phase transi- 
tion from a black hole to a Minkowski embedding takes place. Recall that Ad? ^ NtN^T'^. 




0.2 0.4 0.6 0.8 1 1.2 

m = M/T 



Figure 4: D3/D7 system: Trace of the spectral function as a function of m for (from top to bottom on 
the left-hand side) oj — {0, 0.6, 1.2, 1.4, 1.9, 4, 9}. The top, solid, red curve yields (up to normalisation) 
the electric conductivity. 
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4.2 Massive quarks 

The action (3.6) for tlie D7-branes in the absence of a gauge field takes the form 



(4.22) 



By varying this with respect to tp{u) one obtains a non-Hnear, second-order differential equa- 
tion for the D7-branes embedding. We were unable to solve this equation analytically for 
Mq 7^ 0, but a numerical solution was obtained in [8, 9]. This was done by integrating the 
differential equation from the horizon towards the boundary for different values of ipQ. From 
the value of the solution at the boundary one then reads off m{tpo) and c('i/'o) according to 
(4.6). Using the solution ijj{u) and the ansatz (4.14) in (3.11) a differential equation for F is 
obtained, which can again be integrated numerically. Substituting the result in (4.16) yields 
the trace of the spectral function. Note that in eq. (4.16) one should think of V'o as a function 
of m. 

As shown in [8, 9], thermodynamically stable 
black hole embeddings for the D7-branes exist for 
< ^0 < 0.941, or equivalently for < m < 1.306, 
and in this parameter range m is a monotonically 
increasing function of ^pQ, as exhibited in fig. 5. 
For ipQ > 0.941 or m > 1.306 the free energy of 
the D7-branes is minimised by a Minkowski embed- 
ding, and so a first order phase transition occurs. 
Black hole embeddings with 0.941 < Vo < 0.9621, 
or equivalently with 1.306 < m < 1.3092, are ex- 
pected to be metastable [21, 22], whereas those with 
ipo > 0.9621 have negative specific heat and are 
therefore unstable [9]. 

Converting the differential emission rate (2.3) to the emission rate per unit photon energy 
and using (4.16) gives 



0.4 
0.2 



Figure 5: Dimcnsionlcss ratio m = M/T 
as a function of ipo in the range of values for 
which stable black hole embeddings exist. 



dr 



(1 - ro 



2\3/2 



W 



2-mij 



1)2^^1^(0;, 0)1' 



(4.23) 



Our results for the spectral function and for the photon emission rate for thermodynamically 
stable black hole embeddings are shown in figs. 3, 4, and 6. Some interesting features of these 
plots are as follows. 

The top, solid, red curve in fig. 3 corresponds to the analytic result (4.21) for massless 
quarks. This is a monotonically decreasing function that for large frequencies decays as u~^/^. 
All curves with < m < 0.93 arc also monotonically decreasing functions of the frequency, 
whereas curves with 0.93 ^ m < 1.306 are increasing functions of lu for small lo. At large 
values of to all curves seem to decay with the same power of the frequency as the Mq = curve, 
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Figure 6: D3/D7 system: Photon emission rate as a function of oj for the same values of to as in 
fig. 3, i.e., for (from top to bottom) to = {0, 0.6, 0.85, 0.93, 1.15, 1.25, 1.306}. 

presumably due to the fact that at large lo the scale associated to the quark mass becomes 
irrelevant. The behaviour at small w can also be easily understood. As a; ^ the suppression 
factor (l — V'o)^''^ ™ ^1- (4-16) dominates, and so for small values of uj the spectral density 
decreases monotonically with m, as shown in fig. 4. This effect is also responsible for the 
fact that the height of the peak of the photon emission curves in fig. 6 largely decreases as 
m increases from zero to its maximum value. The suppression factor above has a simple 
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Figure 7: D3/D7 system: Spectral function for non-stable black hole embeddings. The red curve 
with the longest dashes corresponds to "00 = 0.9621, the green curve with intermediate dashes to 
■00 = 0.979, the blue curve with the shortest dashes to ipQ = 0.999996, and the solid, purple curve to 
ipo = 0.9999981. 



geometric origin as the decreasing area of the induced black hole horizon on the D7-branes 
as V'o increases, as can be seen from eq. (4.5). 

Note that the top, solid, red curve in fig. 4, which corresponds to w = 0, gives (up to 
normalisation) the electric conductivity (2.10). Specifically, denoting by h{m) the curve in 
question, one has: 



A{2ttT) duj 



e2 



^=0 4^ 



NfN,Th{m). (4.24) 



Again, the difference between our A^fA^^^, scaling and the A^^ scaling found in [6] reflects the 
difference in the number of electrically charged degrees of freedom. 

At intermediate values of uj the spectral function is not a monotonic function of m, as 
can be seen in fig. 4. In fig. 3 this is reflected in the fact that curves for different values of m 
cross each other around 1 < w < 2. The same behaviour is of course observed in the plot of 
the photon production shown in fig. 6. 

It is also interesting to examine the spectral function for black hole embeddings beyond 
the phase transition, i.e., in the region in which these embeddings are metastable or unstable. 
The results for the spectral function are shown in fig. 7. The most remarkable feature of these 
plots is the appearance of well defined peaks in the spectral function, which become narrower 
and more closely spaced, seemingly approaching delta- functions, as — *■ 1- We will discuss 
the interpretation of this fact in the last section. 
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5. The D4/D6 system 



We now turn to the D4/D6 system, described by the array 

0123456789 
Z:>4 X X X X X (5.1) 
D6 X X X X XXX 

In the decoupHng Umit, the resulting gauge theory is five-dimensional SYM coupled to fun- 
damental hypermultiplets confined to a four-dimensional defect, so we have again d = 3. In 
order to obtain a four-dimensional gauge theory at low energies, one may compactify x^, the 
D4-brane direction orthogonal to the defect, on a circle. If periodic boundary conditions for 
the adjoint fermions are imposed, then super symmetry is preserved and the four-dimensional 
theory thus obtained is non-confining. In this case the appropriate dual gravitational back- 
ground at any temperature is (3.1) with periodically identified. Instead, if antipcriodic 
boundary conditions for the adjoint fermions are imposed, then supersymmetry is broken 
and the four-dimensional theory exhibits confinement [12] and spontaneous chiral symme- 
try breaking [14] . The holographic description at zero-temperature consists then of D6-brane 
probes in a horizon-free background, whose precise form is not needed here. At a temperature 
Tdec set by the radius of compactification, the theory undergoes a first order phase transition 
at which the gluons and the adjoint matter become deconfined. In the dual description the 
low-temperature background is replaced by (3.1). If T^^^ < T(^n, the D6-branes remain outside 
the horizon in a Minkowski embedding. As T is further increased up to Tfu^ a first order phase 
transition for the fundamental matter occurs. 
The metric (3.1) for the black D4-brane is 

ds^ = ^ i-fdxl + cix^) + + r^/^L^/^dn^ , (5.2) 

where 

_ 16 

■> ^3 ' ""o - "9" 

The five-dimensional Yang-Mills coupling constant in the dual gauge theory is dimensionful 
and given by = ^'^'^ds'^s- It is convenient to introduce the dimensionless coordinate 
u = rQ^^/2r^/^, in terms of which the metric becomes 



L^ = TrgsNJl, / = 1 - , ro = -7r^T^L\ (5.3) 



with / = 1 — 4u^. As before, the horizon is at n = 1/2 and the boundary at u ^ 0.^ Since the 
D6-branes wrap a two-sphere in the directions transverse to the D4-branes, it is also useful 

■'The D4-branc metric considered in this section is not asymptoticaUy of the form AdS times a sphere. The 
framework for the calculation of correlators is less well developed for such backgrounds, so we will proceed by 
analogy with the AdS case. Presumably, however, this procedure can be made rigourous by lifting the D4- 
brane geometry to M-thcory, in which it becomes an M5-branc geometry whose asymptotic form is AdSr x S^. 
The D6-brane lifts in turn to a KK-monopole, whose worldvolume effective action can be found in [29]. 



-17- 



to write the metric on the four-sphere as 

dn'] = de^ + sin^ Odnl + cos^ edip^ , (5.5) 

and to set = cos^. The D6-brane embedding is then specified as ^ = il>{u), in terms of 
which the induced metric takes the form 

ds,, = ^-vrTLj - {-fdx, + dx [[ 9^2^(1 _ ^2) ) du + {I - )dn, 

(5.6) 

In eqs. (5.2) and (5.4) x = (x^, . . . , x^), whereas in (5.6) and all equations below in this 
section x = (a;^, . . . , x^). 

As in the previous section, the branes embedding 'ip{u) is determined by extremising the 
action (3.6) with F = 0. Note that in this case there is a non-constant dilaton, 



reflecting the non-conformality of the dual gauge theory. Near the boundary, tp{u) behaves 



as 



ip{u)=m u^/^ + c u^l^ + • • • , (5.8) 



where the dimensionless constants m and c are proportional to the quark mass and condensate, 
respectively [8, 9]: 

^'^^ 2^1 = ^^^^^^^'^"^' ^^-^^ 

05/3 

{O) = -2'>/\HlN,T^,rlc = N,Kg^^{TfT^c, (5.10) 

with 

<7eff(T)2 ^XT = gl^.N^T (5.11) 
the effective 't Hooft coupling at the scale T. In this case, we may write m = M'^/T'^ with 

M^^'(^)'^7.m(^]\ (5.12) 



The scale M is again related to the mass gap in the meson spectrum of the D4/D6 system at 
zero temperature through M ~ 0.233 Mn^os [20]. 

Following the previous section one obtains the equation of motion (3.11) with the func- 
tions P{u), Q{u) given in this case by 

^/(27rr)2 2^1-^2 + 9^2_y.^,2 



^In this limit, the coordinate u used here and the coordinate p used in [8, 9] are related through p^^^ = l/i 
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This equation of motion may be derived from the action (3.12), where 




3r(27r£2)2 



2 = ^iVfTV.T^ , 



(5.14) 



D6 



(f 7r2r2L3) 




(5.15) 



is the normahsation constant introduced in [8, 9]. In this case eq. (4.15) becomes 



hm^ 2Q{u)A* {uj,u)duA{uj,u) = 




(5.16) 



and the trace of the spectral density, eq. (2.8), takes the form 



X';(a;) = -4ImG^(a;) = SX^e (l - ^g) 



2<^ |i^(a;,0)| 



(5.17) 



The action (3.6) for the D6-branes in the absence of a gauge field reduces to 



-^De oc 




(5.18) 



Prom this point one proceeds in complete analogy with the D3/D7 system, i.e., one solves for 
the D6-brane embedding numerically and uses the result to solve, also numerically, for the 
function F{uj,u). Substituting this in (5.17) one obtains the trace of the spectral density. 

In this case the phase transition takes place at ipQ = 0.822 or m = 1.589. Below these 
values black hole embcddings arc stable, whereas for 0.822 < Vo < 0.907, or 1.589 < m < 
1.630, they are presumably metastable. For ipQ > 0.907 black hole embcddings possess a 
negative specific heat and are therefore unstable [9]. The results for stable embeddings are 
plotted in figs. 8, 9 and 10, whereas those for metastable or unstable embeddings are shown in 
fig. 11. Remarkably, these plots share many of the qualitative features of the D3/D7 system. 
We will come back to this in the next section. 

6. Discussion 

Results from heavy ion collision experiments at RHIC indicate that the QGP created there 
behaves as a strongly coupled liquid. In the absence of systematic, non-perturbative methods 
to calculate real-time observables in QCD, it is useful to calculate these observables in gauge 
theories for which a gravity dual is known, the hope being that universal properties of such 
theories may be shared by QCD itself. 

In this paper we have calculated retarded correlators of electromagnetic currents at light- 
like momenta in a class of large- A'^., strongly coupled SU{N^) SYM theories with A^f flavours 
of fundamental matter. From these correlators one may extract several observables associated 
to photon production. Our results are valid to leading order in the electromagnetic coupling 
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Figure 8: D4/D6 system: Trace of the spectral function as a function of uj for (from top 
to bottom on the left-hand side) m — {0,1,1.3,1.45,1.52,1.56,1.589}, or equivalently for ■f/^o — 
{0, 0.449, 0.608, 0.703, 0.755, 0.791, 0.822}. 




0.2 

0.25 0.5 0.75 1 1.25 1.5 

Figure 9: D4/D6 system: Trace of the spectral function as a function of m for (from top to bottom on 
the left-hand side) w — {0, 0.6, 0.97, 1.4, 1.9, 4, 9}. The top, solid, red curve yields (up to normalisation) 
the electric conductivity. 
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UJ = k^/2TTT 




0.86 0.88 0.9 0.92 0.94 0.96 

u; = k°/2TTT 

Figure 10: D4/D6 system: Photon emission rate as a function of uj for tlie same values of m as in 
fig. 8, i.e., for (from top to bottom) m = {0, 1, 1.3, 1.45, 1.52, 1.56, 1.589}. 

constant and in N^/N^, but non-perturbatively in the SU{N^) interactions. At a qualitative 
level, these results exhibit a number of universal features. 

For small values of the quark mass the combination x^fii^)/^ is a featureless, monotoni- 
cally decreasing function of uj, as observed in figs. 3 and 8. Instead, for large enough values 
of m, first increases with oj, then reaches a maximum and finally decays as u!~^ for 
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Figure 11: D4/D6 system: Spectral function for non-stable black hole embeddings. The red curve 
with the longest dashes corresponds to ipo = 0.979, the green curve with intermediate dashes to 
ipo = 0.9999, the blue curve with the shortest dashes to ipo — 0.999996, and the solid, purple curve to 
V-o = 0.9999993. 

large lo. The exponent is however not universal: For the D3/D7 system one has 6 = 1/3, 
whereas for the D4/D6 one finds 6 = 2/3. 

The qualitative behaviour of x%{^)/'-^ as a function of m for fixed values of oj also exhibits 
universal features. For small u) this is a monotonically decreasing function of m, as displayed 
in figs. 4 and 9; in particular, the top, solid, red curves in these figures, which correspond 
to m = 0, yield (up to normalisation) the electric conductivity.'' Instead, for intermediate 
values of oj, the spectral function as a function of m ceases to be monotonic. In figs. 3 and 8 
this is reflected in the fact that curves for different values of m cross each other. This is also 
the case for the photon production curves, as can be seen in figs. 6 and 10. This is perhaps a 
counterintuitive result, as it means that, at certain frequencies, plasmas with heavier quarks 
glow more brightly than those with lighter quarks. 

Above we have calculated the rate of photon emission of an homogeneous, infinitely 
extended, thermally equilibrated plasma. In order to extract a prediction for a physical situ- 
ation such as that of a heavy ion collision experiment, our results would have to be integrated 
over the spacetime evolution of the plasma created in such experiments, as determined by 
some hydrodynamic model [30]. In such a situation the temperature changes over time (and 
space) from some initial, maximum temperature T^^^, to some final, minimum temperature 
Tnjin, and one measures the total number of produced photons of a given energy. This means 

^In the D4/D6 system the electric conductivity is a = {e^ /i-K)N{NcT h(rn), where h{m) is the top, solid, 
red curve in fig. 9. 
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that, unlike in the plots presented above, both u and m vary simultaneously as a; oc and 
m oc T"^ (for the D3/D7 system) or m oc T"^ (for the D4/D6 system). Temperatures around 
r^iax — 4rdee are expected to be achieved at LHC, whereas those at RHIC are roughly around 
2Tdge, where T^ec — 175 MeV is the deconfinement temperature in QCD. On the other hand, 
Tdec is the lowest temperature for which a plasma description is appropriate, since below such 
a temperature hadronisation occurs. For illustrative purposes, let us therefore consider the 
photon emission predicted by our models as a function of T, for Tjec <T< AT^^^ and for 
fixed values of and Mq. It is then convenient to write 



T 

dec -l\ 

uj = -fix, 71 = ——, x = — . (6.1) 

^TTi dsn -£ 



Additionally, for the D3/D7 system we have 



m = j2X, 72 = ^ , (6.2) 



whereas for the D4/D6 system 



dec 

Recall that M^^^ oc is the typical mesonic scale in the gauge theory at zero temperature. 

We are interested in the rate of photon production as x varies between 1/4 and 1 for 
fixed values of 71,72- Converting the differential emission rate (2.3) to the emission rate per 
unit photon energy and using the results (4.16), (5.17) for the spectral functions gives 

a^MNtN^ik^fT^,, dk^ X (62^^71^ - 1) 271^ \F{jix,0)\'^ ' ^ ' ' 

where = 1,^2 = 3/2} for the D3/D7 system, = 3/4,(32 = 1} for the D4/D6 system, 
and the factor (1 — iPq) must be understood as a function of m{x). 

The right-hand side of cq. (6.4) for the D3/D7 system is plotted in fig. 12 for several values 
of the photon energy and the mesonic scale M^^^. Specifically, the top graph corresponds to 
k^ = 100 MeV, whereas the bottom graph corresponds to = 1000 MeV. The meson masses 
correspond to representative values in QCD for bound states of u and d quarks, s quarks 
and c quarks: M^,, = M^(140 MeV), (1020 MeV), Mj/^ (3096 MeV). For comparison, we 
have also included curves for M^^^ = and M^^^ = M* , where M* ~ 0.766Tdcc — 1435 
MeV is the critical mass (corresponding to m = 1.306) for which the phase transition to a 
Minkowski embedding takes place exactly at T = T^^^. For M^^^ = Mjj^ this transition 
takes place at Tf^^ c± 377 MeV > T^^c-, and this is the reason why the dotted, orange curve 
terminates at a; = 175/377 ~ 0.464: At this temperature the corresponding quarks form 
mesonic bound states, i.e., in the string description the probe D7-brancs jump from a black 
hole to a Minkowski embedding. Fig. 13 displays an analogous plot for the D4/D6 system, 
for = 100 MeV. In this case the dashed, green line also terminates at T ~ 188 MeV > T^eci 
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Figure 12: D3/D7 system: Photon emission rate as a function of T for fixed /c" and Mmcs- The top 
graph corresponds to fc° = 100 MeV, whereas the bottom one corresponds to = 1000 MeV. From 
top to bottom, the different curves correspond to Mmcs — 0, 140, 1020, 1435, 3096 MeV. Note that the 
first two curves (solid red and blue dashed) are virtually coincident. 



since for the D4/D6 system both the J ftp and the (p niesons survive as bound states above 
the deconfinement temperature [9]. 

One common feature of both models is that the rate of photon production is most sensitive 
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Xicc / T 

Figure 13: D4/D6 system: Photon emission rate as a function of T for = 100 MeV and several 
values of M^cs: From top to bottom, the different curves correspond to M^es = 0, 140, 947, 1020, 3096 
MeV. Note that the first two curves (solid red and blue dashed) are virtually coincident. 

to the quark mass for low-energy photons and low temperatures. As we see in fig. 12, the 
difference between the curves is greater for = 100 MeV than for = 1000 MeV, and in 
the former case it is greater for low temperatures. We have verified that as increases above 
1000 MeV, the difference between curves becomes smaller and smaller. The same conclusion 
holds for the D4/D6 system. It is also remarkable that, contrary to what one might have 
perhaps expected, the rate of photon production associated to heavier quarks is smaller than 
that of lighter quarks, but it is not negligibly smaller. For example, we see in the figs, above 
that the number of photons emitted by c quarks is a significant fraction of the number of 
photons emitted by u and d quarks for a wide range of temperatures. 

Before leaving our discussion of possible implications for heavy ion collision experiments, 
we would like to emphasise that our results apply to photons produced by a thermally equili- 
brated plasma. However, in heavy ion collisions there are other sources of photon production. 
For example, prompt photons are produced by initial collisions of the partons that constitute 
the colliding ions, whereas decay photons are produced by the decay of certain hadrons once 
the plasma ball has hadronised. Thus comparison of any theoretical prediction for thermal 
photons to empirical data requires being able to distinguish experimentally between these 
different sources. 

The features described above refer to thermodynamically stable black hole embeddings, 
but we also obtained results for metastable and unstable embeddings. These embeddings are 
skipped over by the meson-melting phase transition, but they are still interesting since they 
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illustrate how the spectral function approaches that of Minkowski embeddings, as discussed 
in [22]. Results for the spectral function for several near-critical embeddings, i.e., for values 
of -00 close to 1, are displayed in figs. 7 and 11. Clearly the spectral function develops closely 
spaced, narrow peaks that seemingly approach delta-functions. Thus the form of the spectral 
function appears to approach the form we expect for Minkowski embeddings,^ namely an 
infinite sum of delta functions supported at a discrete set of energies = |k|. Each of these 
delta-functions is associated to a meson mode on the Dq-branes with null momentum. The 
existence of these modes may seem surprising in view of the fact that, as expected on general 
grounds and as verified in ref. [9], the meson spectrum in the Minkowski phase possesses 
a mass gap, but in fact it follows from the analysis in ref. [9]. To sec this, consider the 
dispersion relation k^{U) for a given meson in the Minkowski phase. The fact that there is 
a mass gap means that fe'^ > at k = 0. On the other hand, in the limit of infinite spatial 
momentum, |k| oo, the dispersion relation takes the form ~ v\\l\ with v < 1. The 
reason for this is easily understood: For larger and larger spatial momenta, the wave function 
of the meson becomes more and more concentrated at the tip of the Dq-branes, and so the 
speed of the meson is simply the local speed of light at this lowest point. Because of the 
gravitational redshift, this speed is always subluminal. Continuity then implies that there 
must exist a value of k such that fc'^dkl) = |k|. This is illustrated by the fact that the solid 
black line in fig. 14 intersects the other two curves. Although the curves shown in the fig. 14 
correspond to scalar mesons, it is clear from the arguments above that an analogous result 
would hold for other types of mesons, in particular for vector-like mesons. Since these mesons 
are absolutely stable in the large- iV^, strong coupling limit under consideration, they give 
rise to delta-function- like (i.e., zero- width) peaks in the spectral function of electromagnetic 
currents. 

In this paper we have studied photon production by calculating the electromagnetic 
current-current correlator (2.5) at light-like momenta. It would be interesting to extend the 
calculation to time-like momenta, since this determines the nuumber of dileptons produced by 
virtual photon decay. A first step in this direction was given in ref. [22] , where the correlator 
was calculated for vanishing spatial momentum. 

A further extension consists of calculating the rate of photon and dilepton production in 
the presence of a finite baryon chemical potential [31, 32] of density [28, 33]. For non-zero 
baryon density black hole embeddings exist for all values of T and Mq, and so for large MJT 
the spectral function of electromagnetic currents should again reveal high and narrow peaks 
corresponding to the existence in the spectrum of long-lived vector mesons [32] . 
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Figure 14: Dispersion relation ijj(|fc|) for the lightest spin-zero mesons on a D7-brane Minkowski 
embedding in a D3-brane background [9] . The solid blue curve corresponds to a pseudo-scalar meson, 
whereas the red dashed curve corresponds to a scalar meson. The solid black line corresponds to 
a; = |k|. 
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